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We investigate the hydrodynamic effects on the dynamics of critical concentration fluctuations 
in multicomponent fluid membranes. Two geometrical cases are considered; (i) confined membrane 
case and (ii) supported membrane case. We numerically calculate the wavenumber dependence of 
the effective diffusion coefficient by changing the temperature and/or the thickness of the bulk fluid. 
For some limiting cases, the result is compared with the previously obtained analytical expression. 
An analogy of the multicomponent membrane to 2D microemulsion is explored for the confined 
membrane geometry. 



I. INTRODUCTION 

Biomembranes can be regarded as multipurpose en- 
velopes fundamental to the very existence of life. Com- 
posed of a huge variety of amphiphilic molecules, this ap- 
proximately 5 nm thick quasi-two-dimensional fluid sheet 
separates the inner and outer environments of the cells 
and organelles [1]. Apart from delineating the cell and 
organelle boundaries, membranes also play a significant 
role in a variety of physiological functions such as trans- 
membrane transport of materials and cell signaling 

Biomembranes have attracted renewed interests in the 
context of the lipid "raft" hypothesis proposed a little 
over a decade ago Q. These 10-100 nm sized lipid do- 
mains with higher concentrations of cholesterol were pro- 
posed to play a significant role in regulating certain cel- 
lular functions |2|-|4j ■ Recent experiments using the stim- 
ulated emission depletion (STED) microscopy on plasma 
membranes in vivo narrowed the size-range of the rafts 
to 10-20 nm [f|. Despite extensive studies in this area, 
the details of the underlying physical mechanisms leading 
to formation of rafts, their stability, and the regulation 
of the finite domain size remain elusive and controver- 
sial 0,0]. 

Although the biological significance of rafts is still un- 
der debate, there is no doubt about the rich physics and 
chemistry that have been uncovered by the membrane 
studies to examine the raft hypothesis. Numerous exper- 
iments on intact cells and artificial membranes containing 
saturated lipids, unsaturated lipids and cholesterol have 
demonstrated the segregation of lipids into liquid-ordered 
(L Q ) and liquid-disordered (Lh) phases below the miscibil- 
ity transition temperature [g, |9[ . The L Q -phase is usually 



rich in saturated lipids and cholesterol. Below the transi- 
tion temperature, the domains undergo coarsening with 
the largest domain limited by the system size |9|. The 
primary driving force for the domain coarsening is due to 
the positive line tension at the domain boundaries. Stud- 
ies on the diffusion of domains (Tol j and on the dynamics 
of domain coarsening [iH - flli ] have been reported. 

Recently, studies on multicomponent membranes above 
the transition temperature have also gained much atten- 
tion. As the critical point is approached from above, 
one observes composition fluctuations spanning a wide 
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range of length and time scales [14j. Veatch et al. 
made a notable attempt to investigate critical fluctu- 
ations in lipid mixtures [T3|. Deuterium NMR ex- 
periments on model ternary membranes composed of 
dioleoylphosphatidylcholine (DOPC), dipalmitoylphos- 
phatidylcholine (DPPC) and cholesterol were used to 
construct the ternary phase diagram. With the deter- 
mination of the line of miscibility critical points, they 
observed that the NMR resonances were broadened in 
the vicinity of the critical points. Such a spectral broad- 
ening was attributed to the compositional fluctuations in 
the membrane having spatial dimensions less than 50 nm. 

A more quantitative analysis of the critical fluctua- 
tions using fluorescence microscopy was addressed by 
Honerkamp-Smith et al. for ternary mixtures of DPPC, 
diphytanoylphosphatidylcholine (diPhyPC) and choles- 
terol [15]. When the critical temperature T c is ap- 
proached from above, the correlation length diverges ac- 
cording to £ w \T — T c \~ v where v is the critical expo- 
nent and T the temperature. (In order to prevent the 
confusion with the notations used later, the critical ex- 
ponents are written with a bar.) When T c is approached 
from below, on the other hand, the order parameter 
given by the difference in lipid compositions vanishes as 
S^p ~ (T c — T)P . From the measurements of the critical 



exponents v and /3, the authors concluded that the criti- 
cal behavior in ternary membranes is in the universality 
class of the 2D Ising model [16j . Later, it was also shown 
that giant plasma membrane vesicles extracted from that 
of living rat basophil leukemia cells exhibit a critical 
behavior [13]. These experimental observations suggest 
that lateral heterogeneity present in real cell membranes 
at physiological conditions correspond to critical fluctua- 
tions [111 [IB- In other words, concentration fluctuations 
above the transition temperature (rather than below) can 
also be responsible for raft structures in cell membranes. 

There have also been several theoretical works on con- 
centration fluctuations in multicomponent membranes. 
Using renormalization group techniques, Tserkovnyak 
and Nelson calculated protein diffusion in a multicom- 
ponent membrane close to a rigid substrate [Hj]. They 
pointed out that, in the vicinity of the critical point, the 
effective protein diffusion coefficient acquires a power- 
law behavior. Scki et al. reported equivalent results 
with the use of a two-dimensional (2D) hydrodynamic 
model invol ving a phenomenological momentum decay 
mechanism [19]] • Later Haataja showed that the effec- 
tive diffusion coefficient exhibits a crossover from a log- 
arithmic behavior to an algebraic dependence for larger 
length scales [20j . In his theory, an approximate empiri- 
cal relation for the diffusion coefficient of a moving object 
was employed [2l| . A rigorous hydrodynamic calculation 
performed by Inaura and Pujitani arrived at similar re- 
sults (12]. 

The present article uses the idea of critical phenomena 
to calculate the effective diffusion coefficient in multicom- 
ponent lipid membranes. Based on the time-dependent 
Ginzburg-Landau approach with full hydrodynamics, we 
calculate in particular the decay rate of the concentra- 
tion fluctuations occurring in membranes. We deal with 
the case where the membrane is surrounded by a bulk 
solvent and two walls as depicted in Fig. Q] Such a situ- 
ation is worth considering because biological membranes 
interact strongly with other cells, substrates or even the 
underlying cytoskeleton which affects the structural and 
transport properties of the membrane [23[ . We note that 
effects of the substrates on the diffusion coefficient of pro- 
teins have been investigated before [24|, [25[. Two sur- 
rounding geometries of the membrane are discussed; (i) 
confined membrane and (ii) supported membrane. We 
also study the situation when the multicomponent mem- 
branes form 2D microemulsions [26]. This interesting 
viewpoint is motivated by a recent work which predicts 
the reduction of the line tension in membranes containing 
saturated, unsaturated and hybrid lipids (one tail satu- 
rated and the other unsaturated) [13, HI]. Based on chain 
entropy arguments, they proposed that hybrid lipids in 
2D play an equivalent role to surfactant molecules in 
three-dimensions (3D). Hence we shall explore the con- 
centration fluctuations in 2D microemulsion. 

This paper is organized as follows. In Section [Til we 
first obtain the membrane mobility tensors that are used 
in the subsequent calculations. In Section ITO1 the decay 
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FIG. 1. Schematic picture showing a planar liquid membrane 
having 2D viscosity rj located at z = 0. It is sandwiched by 
a solvent of 3D viscosity 77*. Two impenetrable walls are 
located at z = ±fa bounding the solvent. 



rates of concentration fluctuations are analyzed for con- 
fined and supported membrane cases. Section llVI deals 
with the microemulsion picture of the multicomponent 
membranes. Finally we close with some discussions in 
Section [Vj 



II. MEMBRANE HYDRODYNAMICS 

In this section, we start with the membrane equations 
of motion as well as the requisite membrane mobility ten- 
sors for the two geometries. The membrane is assumed to 
be an infinite planar sheet of liquid and its out-of-plane 
fluctuations are totally neglected, which is justified for 
typical bending rigidities of bilayers. The liquid mem- 
brane, fixed in the rry-plane at z — 0, is embedded in a 
bulk fluid such as water or solvent that is further bounded 
by hard walls as shown in Fig. [T] The upper (z > 0) and 
the lower (z < 0) fluid regions are denoted by " + " and 
" — " , respectively. Since the 3D viscosities of the upper 
and the lower solvent can be different, we denote them as 
r)g . Consider the situation in which impenetrable walls 
are located at z = ±fo , where h + and h~ can also be 
different in general. 

Let v(r) be the 2D velocity of the membrane fluid. 
The 2D vector r = (x, y) represents a point in the plane 
of the membrane which is assumed to be incompressible 

V-v = 0. (1) 

Here V is a 2D differential operator. We work in the 
low-Reynolds number regime of the membrane hydrody- 
namics so that the inertial effects can be neglected. This 
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allows us to use the 2D Stokes equation given by 
??V 2 v - Vp + f s + F = 0, 



(2) 



where rj is the 2D membrane viscosity, p the 2D in-plane 
pressure, f s the force exerted on the membrane by the 
surrounding fluid ("s" stands for the solvent), and F is 
any other force acting on the membrane. 

Once we know f s , the membrane velocity can be ob- 
tained from eqn ([2]) as 



v[q] = G[q]-F[q] 



(3) 



where v[q], G[q] and F[q] (q = (q x ,Qy)) are the Fourier 
components of v(r), G(r) and F(r) defined by 



and 



v ( r ) = J V W ex p(*q • r )> 

G ( r )=y ^^G[q]exp(iq- r), 
F(r) = J ^F[q]exp(iq-r), 



(4) 



(5) 



(6) 



respectively. Using the stick boundary conditions at 
z = and z = , we can solve the hydrodynamic 
equations to obtain f s [2{| . Its perpendicular component 
with respect to q is given by 

/s±[q] = -r]+v ± qcoth(qh + ) - r]~v±qcoth(qh~), (7) 

where q = |q| and v± is also the perpendicular compo- 
nent of v with respect to q. On the other hand, one 
can show that the parallel component of f s is zero. After 
some calculations, the components of the mobility tensor 
G[q] finally becomes 



GW[q] 



l 



T]q 2 + q[rj£ coth(qh+) + r] s coih(qh )] 
q*qp 



x I <W q2 



(8) 



with a, ft = x, y. The full details of the calculation arc 
given in the separate article by Ramachandran et al. [2!| . 

As in ref . [22| , we first consider the case when the two 
walls are located at equal distances from the membrane, 
i.e., h + = h~ = h. Then the above mobility tensor 
simplifies to 



•q[q 2 + vq coth(g/i)] 



>a/3 



q a qi3 



(9) 



where v = 2rj s /rj with r/ s = (77+ +rj~)/2. This expression, 
taken in the limit of large h, has been used previously in 
calculating the correlated diffusion coefficients of parti- 
cles embedded in a membrane [301 ] . 



For supported membranes, h + is infinitely large while 
h~ is finite when compared to the membrane thickness. 
In this case, eqn (jHJ reduces to 29] 



G a p [q] 



f 



?7[g 2 -(-^(l + coth( g /i-))/2] 



q a qfi 

q 2 



(10) 

This equation has been used for the investigation of the 
correlated dynamics of inclusions in a supported mem- 
brane 1311] . 



III. DYNAMICS OF CONCENTRATION 
FLUCTUATIONS 

In order to discuss the dynamics of concentration fluc- 
tuations above the transition temperature, we closely fol- 
low the formalism used in ref. [32j | . Here we extend our 
previous work in ref. [l9| for membranes embedded in 
a solvent confined by two walls (confined membrane) or 
supported on a substrate (supported membrane). 

Consider a two-component fluid membrane composed 
of lipid A and lipid B whose local area fractions are de- 
noted by </>a(f) and 0b (r), respectively. Since the rela- 
tion 4>a(t)+4>b(^) = 1 holds, we introduce a new variable 
defined by i/j(r) = 0a (r) — 0b (r). The simplest form of 
the free energy functional J~{ip} describing the fluctua- 
tion around the homogeneous state is 



= / dr 



^ 2 + ^(v^) 2 



(11) 



where a > is proportional to the temperature difference 
from the critical temperature, c > is related to the line 
tension and \x is the chemical potential. 

The time evolution of concentration in the presence 
of hydrodynamic flow is given by the time-dependent 
Ginzburg-Landau equation for a conserved order param- 
eter [11 



dip 
~dt 



V • (yip) = LV 



,5T 



(12) 



where L is the kinetic coefficient. In the membrane hy- 
drodynamic equation @, on the other hand, we need to 
incorporate the thermodynamic force due to the concen- 
tration fluctuations. Hence we have 



-V>v 



ST 

1$' 



(13) 



We implicitly assumed that the relaxation of the ve- 
locity v is much faster than that of concentration ip [r| . 
The membrane velocity can be formally solved using the 
appropriate 2D mobility tensor G a p(r,r') derived in the 
previous section, 



v a (r,t) = \ dr'G^(r,r')(V^) 



ST 
5ip(r') 



(14) 
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Since our interest is in the concentration fluctuations 
around the homogeneous state, we define 6ip(r,t) — 
tp(r,t) — ip, where the bar indicates the spatial average. 
The free energy functional expanded in powers of Sip be- 
comes, 



= j dr[|(5V) 2 + |W) 2 



(15) 



Substituting eqn (|T4|) into eqn (fT2|) . we get 
ST 



d5ip{v,t) = 2 



dt S(Sip) 
dr'(V Q ^(r))G^(r,r')(V^^(r'))^^y. 



ST 



(16) 



We now consider the dynamics of the time-correlation 
function defined by 



S(v,t) = (6^(r 1 ,t)Sip(r 2 ,0)}, 



(17) 



where r = T2 — ri. Within the factorization approxima- 
tion [32| . the spatial Fourier transform of S(r,t) defined 
by 



S(r,t) 



dq 



(2tt)2 

satisfies the following equation 
dS[q,t] 



S[q,t]exp(iq- r), 



dt 



(r«[q]+r( 2 )[q]) S[q,t]. 



(18) 



(19) 



In the above, the first term rW[q] denotes the van Hove 
part of the relaxation rate given by 

r( 1 )[q] = Lfc B Tg 2 X - 1 [q]. (20) 
Here the static correlation function is defined by 

k B T 



X[q] = (<5V[q]^[-q]> = 



c(g 2 +C- 2 )' 



(21) 



where £ = (c/a) 1 / 2 is the correlation length, k-Q the Boltz- 
mann constant, and T the temperature. 

As for the second term in eqn (fl"9|. r( 2 '[q] denotes the 
hydrodynamic part of the decay rate 

r(2)[q] = xR / 7CT 9 « G ^W«^ +P ^ (22) 

where either eqn © or eqn (fTUl) will be used for the 
mobility tensor G a p. 



A. Confined membrane 

First we consider the situation in which the membrane 
is confined by two walls which are located at equidis- 
tant h from the membrane. The appropriate form of the 
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FIG. 2. Scaled effective diffusion coefficient D as a function of 
Q for H = 0.01, 1, 100 when X — 1 for the confined membrane 
case. The numerically calculated points are joined by lines for 
clarity. The solid red lines are from the analytical expression 
given in eqn (|26|) obtained in the limit of small H. 



mobility tensor is given by eqn ([9]). Substituting it into 
eqn (|22p , the hydrodynamic part of the decay rate is ex- 
pressed as 



vxm J 



dp 

(2tt)2 

x\p] 



2 2 



|p - q| coth(|p - q\h) 



|p-q|- 



(23) 



We introduce an effective diffusion coefficient D[q\ (due 
only to the hydrodynamic effect) defined by 



r( 2 )[q]= g 2 D[q] 



(24) 



In order to deal with dimensionless quantities, we rescale 
all the lengths by the hydrodynamic screening length 
v~ l = rj/(2rj s ) such that P = p/v, Q = q/v, X = £v 
and H = hv. Then D[q\ can be rewritten as 

D[Q;X,H] = ^f(l + Q 2 X 2 ) 



poo p2 

/ dP 

Jo Jo 



dd 



P A sin 2 



(1 + P 2 X 2 )[G 2 + G 3 / 2 coth(VG#)] ' 

(25) 



with G = P 2 + Q 2 — 2PQ cos 9. Since this integral cannot 
be performed analytically, we evaluate it numerically. We 
explore the dependencies of D on the variable Q, and the 
parameters X and H . 
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FIG. 3. Scaled effective diffusion coefficient D as a function of 
Q for X = 0.01, 1, 100 when H = 1 for the confined membrane 
case. 
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FIG. 4. Scaled effective diffusion coefficient D as a function 
of X for # = 0.01, 1, 100 when Q = 10" 3 for the confined 
membrane case. The solid red lines are from the analytical 
expression given in eqn (1261) obtained in the limit of small H . 



In Fig. [2j we plot the diffusion coefficient D (scaled 
by k^T /Ajtri) as a function of dimensionless wavenumber 
Q for different solvent thickness H while the correlation 
length is fixed to X = 1 (i.e., fixed temperature). In the 
limit of Q <C 1, D is almost a constant. The calculated 
D starts to increase around Q « 1 and a logarithmic 
behavior (extracted via numerical fitting) is seen for Q 3> 
1. This trend has been previously reported by Seki et 
al. [l^l or Inaura and Fujitani [22 |. When H is small 
such as H = 0.01, the value of D decreases by one order 
of magnitude compared to H = 100. Figure [3] shows the 
diffusion coefficient I? as a function of wave number Q for 
different X (i.e., different temperature) while the solvent 
height is fixed to H = 1. Again, D is nearly constant for 
Q <C 1, and follows an S'-shaped curve with increasing 
Q. Finally, a logarithmic dependence is observed from 
numerical fitting. We note that the above logarithmic 
behavior for Q ^ 1 is in contrast to that of 3D critical 
fluids as given by the Kawasaki function which increases 
linearly with q [33|. 

In Fig. |U we explore the effect of the correlation length 
X on D for different values of H when Q = 10 -3 . The 
quantity X measures an effective size of the correlated re- 
gion formed transiently in the membrane due to thermal 
fluctuations. When X -C 1, the diffusion coefficient D 
decreases o nly logarithmically, which is typical for a pure 
2D system [34M36I ] . When X 1, on the other hand, the 
behavior of D depends on the value of H . The proxim- 
ity to the walls results in a loss of momentum from the 
membrane [l9|, [29|, [37[ ■ This leads to a rapid suppression 
of the velocity field within the membrane such that the 



concentration fluctuations decay slowly. Consequently, 
the values of D are lower for smaller H . The flattening 
of the curves for large X is due to the dominance of the 
X 2 terms in the numerator and denominator of eqn (|25[) . 
In Fig. [SJ we plot D as a function of H for different values 
of X when Q = 10~ 3 . In general, there is a monotonic 
increase in D with larger H followed by a saturation to 
a constant value. We see that the effect of H is most 
prominent for large X (close to the critical point), while 
there is only a weak dependence for I < 1 (far from 
the critical point). The former reflects the fact that the 
membrane fluid is affected by the the outer environment 
when the correlation length £ is larger than the hydrody- 
namic screening length v~ x i.e, X ^> 1 [29j]. For correla- 
tion lengths smaller than the v ~ , the outer environment 
surrounding the membrane is less important. In this sit- 
uation, the system behaves essentially as a pure 2D one. 
When the correlation length becomes larger than v~ l , as 
is the case near the critical point, the outer environment 
significantly affects membrane dynamics. 

For small h, the vq coth(q/i) term in the mobility ten- 
sor eqn ([9]) can be replaced by a constant vjh. In this 
case, Seki et al. obtained an analytical expression for the 
effective diffusion coefficient 19]. Their result in terms 
of the dimensionless quantities Q, X and H can be re- 
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FIG. 5. Scaled effective diffusion coefficient D as a function 
of H for X = 0.01, 1, 100 when Q = 10" 3 for the confined 
membrane case. 




Q 

FIG. 6. Scaled effective diffusion coefficient D as a function 
of Q for H~ = 0.01, 1, 100 when X = 1 for the supported 
membrane case plotted with filled symbols. For comparison, 
the data in Fig. are also plotted with open symbols. 



produced as 



D[Q;X,H] 



k B T 1 + Q 2 X 2 
Anri 2Q 2 X 2 



In 



X 
H 



+ ^{l + Q 2 X 2 )\n 



where 



and 



H 

X 2 

HQ 

V 2X 2 

Q = 



X 



■In 



H(1 + Q 2 X 2 ) 
Q\ + Q 2 _+Q 2 + n" 

n-Q 2 _-i 



(26) 



(Q 2 X 2 + X 2 /H -if +AQ 2 X 2 , (27) 



Q± = y/Q 2 x 2 ±x 2 /H. 



(28) 



Equation (f2"6")l is plotted using red curves in Fig. [2] for 
H = 0.01 and 1 with X = 1. For H = 0.01, the an- 
alytical and numerical data coincide giving credence to 
accuracy of the numerical solutions. It is seen that even 
for H = 1 the agreement is still acceptable. For H = 100, 
however, a significant deviation is observed (not shown), 
which is expected as this limit is beyond the valid range 
of eqn (|26|). 

The red curves in Fig. @] also represent the analytical 
result of eqn (|2"6"|) . It is seen that the analytical and the 
numerical data points almost coincide for H = 0.01 and 
H = 1. For H = 100, there is significant deviation from 
the numerical data (not shown). However, the agreement 
between the numerical result and the analytical expres- 
sion is beyond the expected range of H <C 1 and reaches 



up to H ~ 1, as pointed out by Stone and Ajdari |25| . 
Hence eqn (|26l) can be useful in analyzing the experimen- 
tal data in many situations. 



B. Supported membrane 

Apart from studying membrane dynamics on vesicles, 
experiments are also conducted on supported membranes 
which have a benefit of avoiding curvature effects [38| . 
We proceed now to calculate concentration fluctuations 
for the supported membrane case. As given by eqn (|10l) . 
the mobility tensor of a supported membrane is slightly 
modified from that of a confined membrane. By substi- 
tuting eqn (|10|) into eqn (|22|) . we obtain the expression 
for the effective diffusion coefficient as 

D[Q- X } H-] = + Q 2 X 2 ) dP J o d9 



P 3 sin 2 6 



(1 + P 2 X 2 )[G 2 + G 3 / 2 (l + coth(VGF-))/2] 



(29) 



where H~ = hrv and G — P 2 + Q 2 — 2PQ cos 9 as before. 

Performing numerical integrations, we plot in Fig. [6] 
the effective diffusion coefficient D as a function of Q for 
different values of H~ when X = 1 (closed symbols). In 
general, the behavior of D is similar to that of the con- 
fined membrane case. When H~ <C 1, the calculated D 
is slightly larger than the confined membrane compared 
with the same value of H . This is because the supported 
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membrane is subjected to only one wall, while the con- 
fined membrane is sandwiched by two walls on both sides. 
The confined and the supported membranes show an al- 
most identical behavior when both H and H ~ are large 
enough, as it should be. The other dependencies of D 
are similar to the confined membrane case except for a 
slight increase at small H~ values. 



IV. MEMBRANE AS A 2D MICROEMULSION 

The role of surfactant molecules in 3D microemulsions 
is to reduce the surface tension at the interface between 
oil and water. In an analogy to 3D microemulsions, hy- 
brid lipids (one chain unsaturated and the other satu- 
rated) act as lineactant molecules which stabilize finite 
sized domains in 2D. In other words, hybrid lipids play 
a similar role to surfactant molecules at the interface be- 
tween L Q and Ld domains. It should be also noticed that 
hybrid lipids form a major percentage of all naturally ex- 
isting lipids [1^, H(| ■ Based on a simple model of hybrid 
lipids, Brewster et at showed that finite sized domains 
can be formed in equilibrium [27], l28| . A subsequent 
model predicted stabilized domains even in a system of 
saturated/hybrid/cholesterol lipid membranes [41(. Be- 
ing motivated by this idea, we calculate the decay rate of 
concentration fluctuations when the free energy of the 
multicomponent membrane has the form of a 2D mi- 
crocmulsion. Here we consider only the confined mem- 
brane geometry, and use eqn ([9]) for the mobility tensor. 

The free energy functional for a microemulsion includes 
a higher order derivative term and is expressed in terms 
of Sip as 



■FmeW} = J dr 



lm 2 



|(v 2 ^) 2 



(30) 

with a, g > and c < 0. The negative value of c creates 
2D interfaces, while the term with positive g is a stabi- 
lizing term. This form of the free energy has been used 
previously to study coupled modulated bilayers [42j]. As 
in the previous section, the decay rate of the correlation 
function can be split into two parts. First, the van Hove 
part becomes now 



Lk B Tq 2 XME [q\ 



(31) 



where L is the kinetic coefficient assumed to be same as 
before, and the static correlation function XME[q] is [H| 



9<t 



By defining 



9o = ~v 



(32) 



(33) 



(34) 
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FIG. 7. Scaled effective diffusion coefficient Dme as a func- 
tion of Q for H = 0.01, 1, 100 when Q = E = 1 for the 
confined membrane case. 



the static correlation function can be also written as 

k B T 



XME[qj 



g[(q 2 -<7o) 2 + 



(35) 



On plotting xme as a function of q, a peak appears at 
q = qo followed by a q~ -decay. The width of the peak 
is given by a, and a lamellar phase appears when a = 0. 
Notice that c = is called the Lifshitz point at which the 
peak occurs for q = [HI]. Using the form of eqn (|35l) . 
we can rewrite eqn (|31[) as 

r«[q]=L OT 2 [(g 2 - 9o 2 ) 2 + a 4 ]. (36) 

As in the previous section, we next write the hydro- 
dynamic part of the decay rate in terms of the effective 
diffusion coefficient DmeM as 

r^[q] =q 2 D ME [q}. (37) 
Using eqn §§§ for the mobility tensor, we obtain Dme as 
knT 



Dme [Q] Qo, H] 



47T 2 ?7 



Jo 



,.2n 

■ j d(9 



[{Q - Qo 



P 3 sin 2 1 



2\2 



f C 

Jo 



dP 



[{P 2 - Q 2 ) 2 + S 4 ][G 2 + G 3 / 2 coth(VGtf)] ' 



(38) 



where P = p/v, Q = q/v, Qo = 1ol v , £ = a l v i H — hv, 
and G = P 2 + Q 2 - 2PQ cos 9. 

In Fig. [71 we plot -Dme as a function of Q for different 
values of H when Qo = ^ = 1 are fixed. When Q <C 
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1, D shows a constant value. We also observe the 2D 
characteristic of logarithmic behavior of D for Q 1. 
For Q <C 1 , the effect of the outer environment is felt with 
the suppression of the diffusion coefficient with smaller 
H. The curves almost overlap when Q > 1 indicating 
the negligible effect of the outer environment at large 
wave numbers. An interesting feature of -Dme is the dip 
occurring at Q sa Q which does not exist for binary 
critical fluids. This can be attributed to the peak at 
q = q in Xme[q] [IE IH; EH- F° r 3D microemulsions, 
however, it is known that the effective diffusion coefficient 
varies linearly with q for large wave numbers [32| . 

Although the analogy between a 3D and 2D mi- 
croemulsion has been invoked, it should be pointed 
out that a 3D microemulsion formed from an 
oil/water/surfactant mixture arises predominantly due 
to the differences in relative affinity between the com- 
ponents. For the 2D case, it is the physical interactions 
from the hydrocarbon chain packing requirements at the 
L /L{j interface that gives rise to the lineactant proper- 
ties of the hybrid lipid. 



V. DISCUSSION 

In summary, we have calculated the decay rate of con- 
centration fluctuations in a multicomponent fluid mem- 
brane for two geometries. First we considered the mem- 
brane surrounded by solvent of finite depth which is fur- 
ther bounded by two walls. In the second geometry, we 
allow the solvent depth on one side of the membrane to be 
infinitely large. This is equivalent to a supported mem- 
brane in experimental situations. The resulting integrals 
for the effective diffusion coefficient of the concentration 
fluctuations are calculated numerically to determine the 
various dependencies. We have also explored a possibil- 
ity of considering the multicomponent membrane as a 2D 
microemulsion. 

The present work follows Seki et al. who were able 
to use analytical means to obtain the effective diffusion 
coefficient [19| . Although their result should be valid only 
in the limit of very small h, the agreement between the 
numerical result and the analytical expression is fairly 
well as long as h < v^ 1 . The opposite limit of very large 
h was studied by Inaura and Fujitani through numerical 
methods (22[. Using the general mobility tensor given 
by eqn (|5J|, we are able to probe all the intermediate 
situations of finite h. We have verified that our results 
properly interpolates between these previous works in the 
limits of H and H —> oo. 

In our case, the measure of the size of the transient 
structures is given by the correlation length £. Close to 
the critical point, £ diverges and the hydrodynamic ef- 
fects of the outer fluid play a significant role in altering 
the diffusion coefficient. It has been previously shown 
through explicit calculations that the diffusion coefficient 
has a logarithmic behavior D f=s ln(l/i?) when the size 
of the diffusing object R is less than the hydrodynamic 



screening length, i.e., R < v~ x [HHl]. When R > v~ x , 
on the other hand, the fluid flow in the bulk leads to the 
diffusion coefficient to show l/i?-behavior (this condition 
occurs when there are no walls) (36j . In the presence 
of w alls o r a substrate, the screening length is altered 
to yjh/u, where h is the distance of the walls from the 
membrane |25|. When R ^> y/h/v, the diffusion coeffi- 
cient now shows an algebraic decay D ss 1/R 2 119]. This 
change is attributed to the loss of momentum from the 
membrane to the walls whereas momentum is conserved 
otherwise (37} . 

In the present study, we have kept the model as simple 
as possible. The bending stiffness of typical membranes is 
of the order of 10 fee T which is sufficiently large enough to 
neglect the out-of-plane displacements of the membrane 
itself. The dynamics of the out-of-plane fluctuations have 
been previously studied by Levine and MacKintosh (4|| . 
In our treatment, we have also neglected the effects of 
membrane curvature which can be significant when the 
radius of curvature becomes close to the hydrodynamic 
screening length. Recent calculations by Henle et al. con- 
sidered the diffusion of a point object on a spherically 
closed membrane [U EH . The extension of this work to 
finite sized objects is a particularly difficult proposition. 

From the experiments on model multicomponent vesi- 
cles, the static critical exponents for the order parameter 
and correlation length were found to have values close 
to j3 = 1/8 and v = 1, respectively Furthermore, 
experiments on giant plasma membrane vesicles mea- 
sured the critical exponent 7 = 7/4 which characterizes 
the critical behavior of the osmotic compressibility fl7j . 
These static exponents coincide with the exact results of 
the 2D Ising model [50| • The description presented 
in this paper uses the mean-field approach and therefore 
the corresponding static exponents are f3 = 1/2, v = 1/2, 
and 7=1, respectively. However, attributing biological 
relevance to these critical fluctuations should carried out 
with caution. This is because the plasma membranes 
do not show phase separation phenomena without chem- 
ical treatment 51]. Extraction of membranes from real 
cells also ruptures the association with the underlying cy- 
toskeleton and other active cellular processes like vesicu- 
lar trafficking. 

As mentioned earlier, experiments on real plasma 
membranes have suggested that the cell maintains the 
membranes at a critical composition [l7|. This leads to 
nanometer-sized composition fluctuations at the physio- 
logical temperatures, although these structures are much 
smaller than what can be resolved through optical mi- 
croscopy. We therefore speculate that there is some 
biological relevance in studying concentration fluctua- 
tions. However, this article has mainly concentrated on 
the dynamics towards the equilibrium state of lipid bi- 
layer membranes. In real cells, there are many active 
non-equilibrium cellular processes that are involved in 
the proper functioning of the cells. It has indeed been 
proposed that nano-domain formation may be related 
to the underlying cytoskeleton [52[. There have been 
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several other models which make use of the active non- 
equilibrium phenomena to explain the existence of finite 
sized domains in multicomponent membranes. A review 
of the non-equilibrium models can be found in the review 
article [Hj]. 

Another experimental system which can be used to 
verify our model is the Langmuir monolayer setup. The 
upper "+" region in a Langmuir monolayer system is 
occupied by air and hence 77+ should be set to zero in 
cqn ||5J). In this case, the mobility tensor is given by 
eqn @, but the definition should be replaced with v — 
T)~ /r). The diffusion coefficient as a function of wave 



vector can be obtained via light scattering techniques. 
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